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CO ■ ABSTRACT 

The detections of both X-ray and radio emission from the cluster Gl in M31 have provided strong 
support for existing dynamical evidence for an intermediate mass black hole (IMBH) of mass 1.8 ± 
0.5 x 1O 4 M at the cluster center. However, given the relatively low significance and astrometric 
accuracy of the radio detection, and the non-simultaneity of the X-ray and radio measurements, this 
identification required further confirmation. Here we present deep, high angular resolution, strictly 



simultaneous X-ray and radio observations of Gl. While the X-ray emission (Lx = 1. 74^.44 x 
10 36 (d/750 kpc) 2 erg s _1 in the 0.5-10 keV band) remained fully consistent with previous observations, 
we detected no radio emission from the cluster center down to a 3<r upper limit of 4. 7 //Jy beam -1 . 
Our favored explanation for the previous radio detection is flaring activity from a black hole low 
mass X-ray binary (LMXB). We performed a new regression of the Fundamental Plane of black hole 
activity, valid for determining black hole mass from radio and X-ray observations of sub-Eddington 
black holes, finding logM BH = (1.638 ± 0.070) logL R - (1.136 ± 0.077) logL x - (6.863 ± 0.790), with 
an empirically-determined uncertainty of 0.44 dex. This constrains the mass of the X-ray source in 
Gl, if a black hole, to be < 9.7 x 1O 3 M at 95% confidence, suggesting that it is a persistent LMXB. 
This annuls what was previously the most convincing evidence from radiation for an IMBH in the 
Local Group, though the evidence for an IMBH in Gl from velocity dispersion measurements remains 
unaffected by these results. 

Subject headings: Accretion, accretion disks — globular clusters: individual: (Mayall-II = Gl) — 
X-rays: binaries — Black hole physics — radio continuum: general 
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1. INTRODUCTION 

The existence of intermediate-mass black holes 
(IMBHs), while theor etically well motivated (e.g. 
IVolonteri fc Per na] 120051 ), remains to be definitively con- 
firmed. With masses 50M© < M B h < 1O 6 M , such 
objects would bridge the gap between the dynamically- 
confirmed stellar mass black holes in the Galaxy 
(McClintock & Remillard 2006) and t he supermassive 
black holes at the centers of galaxies ([Magorrian et al.l 
[19981) . 

IMBHs are believed to form via the deaths o f 
primordial (Population III) stars ([Fryer et al.l l2001f h 
or by direct collapse of la rge amounts of gas (e.g. 
IVolonteri fc Begelmanl I2010D . Alternatively, IMBHs 
can form in the cores of dense star clusters, either 
through successive merge rs of stellar-mass black holes 
(Mill er fc Hamilto n 2002), or via a runaway collision of 
massi ve stars in the cluster center followed by direct co l- 
lapse ([Portegies Zwart et al.l l2004a; Gi irkan et al.l l2004). 
Thus, dense star clusters should be among the best lo- 
cations to search for IMBHs, either via their effects on 
the dynamics of the surrounding stars, or via radiative 
signatures of the accretion of gas onto the black hole. 

Searching for IMBHs using dynamical studies of star 
clusters is only possible within the Local Group (e.g. 
iSeth et aTl 120101 ), where we have sufficient angular res- 
olution to resolve the black hole sphere of influence. 
Arguably the best IMBH candidate identified to date 
via stellar dynamics is at the center of the cluster Gl 
Mayall-II) in M31 . With a mass of (7-17)xlO 6 M 
Mevlan et al.ll200l| ), Gl is one of the most massive star 
clusters in the Local Group, and has been suggested to 
be the tidally s tripped core of a former dwarf galaxy 
([Ma et al.ll2007[ ). Fro m modeling the stellar p hotometry 
and kinematics of Gl, IGebhardt et al.l ([2002f ) suggested 
the presence of a central IMBH of mass 2.01q § x 10 4 M o . 
This claim was challenged by Baumg ardt et al.l ([2003bf ). 
who attributed the high central velocity dispersion to a 
central co ncentration of ste l lar re mnants. However, using 
new data, IGebhardt et al.l (|2005l ) ruled out the no-black 
hole model at 97% significance, and refined the mass es- 
timate to 1.8 ± 0.5 x 10 4 M o . 

The presence of an IMBH in Gl would imply 
the possibility of a direct detection of its accre- 
tion signature. Recent XMM- Newton observations 
dTrudolvubov fc Priedhorskvl[2004HPoolev fc Rappaportl 
120061 ) detected X-rav emission from Gl t hat was pos i- 
tionally coincident with the cluster center (Kong 2007), 
a localization subsequently confirmed to within 0.15" by 
higher-precision astrometry with Chandra and the Hub- 
ble Space Telescope (HST) (jKong et al.ll2010l) . The X- 
ray spectrum and luminosity were consistent between 
the XMM-Newton and Chandra observations, and, for 
a dis tance of 750 kpc ([Vilardell et al.l 120101 : iRiess et all 
120121 ) correspond to an unabsorbed 0.3-7 keV luminosity 
of 2.lt°ol x lO^ergs -1 ([Kong et al.ll2010f) By co nsider- 
ing several possible scenarios, iKong et al.l (|201Q[ ) found 
that only an IMBH accreting from a companion star, 
or an ordinary low mass X-ray binary (LMXB) could 
explain the X-ray observations, although neither X-ray 
spectra nor X-ray astrometry could distinguish between 
these two cases. 



From the observed correlation among radio luminos- 
ity, X-ray luminosity and mass of low-power accreting 
black ho les (the Fundamental Plane (FP) of b lack hole 
activity; iMerloni et all 120031 : iFalcke et aLll2004f ). the ra- 
dio emission from an IMBH is expected to be several 
hundred times greater than that from an LMXB, pro- 
vidin g an ideal discriminant bet ween these two scenar- 
ios ([Maccarone fc Kor ding 20Q6|). The detection by the 
Very Large Array (VLA) of a 28 ± 6 /iJy sou rce within 
1.3" of the center of Gl (jUlvestad et aLM2QQ7l hereafter 
UGH07) was consistent with the radio emission expected 
from an accreting IMBH with the observed X-ray lumi- 
nosity. However, the relatively compact VLA configura- 
tion used meant that the astrometric accuracy was insuf- 
ficient to conclusively associate the radio detection with 
the X-ray source. The 1.3" positional offset, while for- 
mally within the 95% confidence X-ray error circle, leaves 
open the possibility that the two sources are unrelated. 
Furthermore, since LMXBs are known to be variable in 
both the radio and X-ray bands, the non-simultaneity 
of existing X-ray and radio data implies that the LMXB 
interpretation cannot be ruled out. To address these con- 
cerns, we report here the results of deep, high-resolution, 
strictly-simultaneous X-ray and radio observations of Gl, 
using Chandra and the Karl G. Jansky Very Large Array 
(J VLA). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Chandra 

Chandra observed Gl on 2011 June 26 (observation 
12291) using the ACIS-S3 detector. As there was no 
evidence for background flaring, we analyzed the full 
34587 s observation using CIAO 4.3 with caldb 4.4.3. Us- 
ing sources detected by the wavelet detection algorithm 
(CIAO wavdetect) and excluding Gl, we registered the 
observation to the absolute astrometri c frame used for a 
previous Chandra localization of Gl ([Kong et al.ll2010f ). 
As Gl was only 0.'26 off- axis for our Chandra observa- 
tion, we analyzed all events in a circle with a radius 50% 
larger than that enclosing 90% of the source's photons, 
and derived the background in a nearby annulus with 5 
times the area of the source aperture. 

Although a variety of spectral models were fit to the 
X-ray emission from Gl we only report the results for 
an absorbed power-law model; a power-law spectrum is 
expected for either an IMBH accreting at a low fraction 
of its Eddington luminosity or a low-luminosity LMXB. 
For the reported best-fit (x 2 = 1.6 for 3 degrees of free- 
dom) , we performed x 2 fitting after applying a minimum 
bin size of 15 counts and Gehrel's weighting; altering the 
bin size or fitting with the Cash statistic did not strongly 
affect the best-fit parameters. We set the absorption col- 
umn density to the Galactic value of 6.5 x 10 20 cm -2 
and applied the Tubingen-Bo ulder absorpt i on (t babs) 
model with abundances from I Wilms et al.l (120001) and 
photo electric absorption cross-sections from lVerner et all 
(fl996h . 

2.2. J VLA 

On 2 011 June 26, we o bserved Gl for 10 h with the 
J VLA (jPerlev et al.ll201lh . giving a total on-source time 
of 446 min. The bandwidth was split into two 1024-MHz 
basebands centered at 5.0 and 7.45 GHz. Each baseband 
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Table 1 

Unabsorbed fluxes and luminosities of Gl for 
different X-ray bands, calculated assuming a 
source distance of 750 kpc. 



Band Unabsorbed flux Luminosity 
(keV) (lCT 14 erg cm" 2 s -1 ) (10 36 ergs" 1 ) 



0.5-10 


2 5Q+ - 79 


i 74 +o.53 

' '^-0.44 


2-10 


7? +0.78 
Lm ' '-0.63 


1 1Q+0.53 
1 - iy -0.42 


0.3-7 


2 sn +0 - 46 


i 55+O.3I 


0.3-10 


9 7Q+0-73 
z - ' y -0.63 


oc+0.49 
J-So-0.43 



consisted of eight 128-MHz sub-bands, each comprising 
64 spectral channels of width 2 MHz. We observed in full 
polarization mode, using an integration time of 3 s. The 
JVLA was in its most-extended A-configuration, giving 
a maximum resolution of 0.23". 

To render our measurement insensitive to artefacts 
gen erated at th e phase center by correlation errors (see, 
e.g. lEkersl[T999h . we pointed 2.5" south of Gl. The am- 
plitude scale was set via observations of 3C48, using the 
coefficients derived at the JVLA by NRAO staff in 2010. 
The secondary calibrator J0038+4137 was observed for 
one minute of every six to permit amplitude and phase 
calibration. Every hour, we observed the nearby check 
source, J0025+3919, to estimate coherence losses via 
comparison of its measured flux density before and af- 
ter self-calibration. 

Data were reduced following s tandard procedures 
within CAS A (Mc Mullin et al.l 120071 ). Imaging was car- 
ried out using natural weighting, for maximum sensi- 
tivity. We deconvolved bright sources out to twice the 
distance of the half-power point of the primary beam 
at 5 GHz, using the ^projection algorithm to prevent 
phase errors due to the sky curvature from affecting our 
deconvolution. 

3. RESULTS 
3.1. Chandra 

With 88.8 net counts, our best fit spectrum for data 
in the 0.3-10 keV range had a power-law index of T = 
1.591q'JJ, where all error bars from X-ray spectral fitting 
indicate 90% confidence intervals. We found a centroid 
position for Gl of 00 h 32 m 46 s .536, 39°34'40".51 (J2000), 
with a la positional uncertainty ellipse having semi- 
major and semi- minor axes of 0.18 and 0.16 arcsec, re- 
spectively. The uncertainties are dominated by system- 
atic errors from the registration of the image to the ab- 
solute astrometric frame. 

We measured an unabsorbed 0.5-10 keV flux of 
2.59^0*66 x 10 -14 ergcm -2 s -1 , corresponding to a lumi- 
nosity of 1.74tg;g x 10 36 (d/750 kpc) 2 ergs" 1 . For com- 
parison, we provide in Table Q] the unabsorbed fluxes and 
luminosities for the different energy bands used in the 
literature both for Gl and to derive the FP regression 
coefficients. 

3.2. JVLA 

No radio source was significantly detected within 28" 
of the center of Gl. Our final image (Figure [1]), 
made using the data from both basebands, had a noise 




0Q h 32 m 47 fi .2 46 3 .S 46*. 5 46\+ 46 s . 2 46 S .Q 
J20O0 Right Ascension 



Figure 1. Radio images of the field containing Gl. Greyscale 
shows the JVLA image from 2011 June 26. Contours show our 
images of the archival VLA data (black contours for 2006 November 
24-26, white contours for 2007 January 14-15), with contours at 
intervals of ±2,4 times the rms noise (7.2 and 16.3 fiJy beam -1 , 
respectively). Blue and cyan error circles denote the core and half- 
mass r adius of the cluster, respectively (position f rom t he HST 
data of IKong et alJ ([20101 ) and radii from Ma et al. (2007)), gree n 
circle denotes the 95% confidence error circle of Kong et al. ( 2010), 
and red ellipse denotes the 95% confidence error ellipse from our 
Chandra observations. Beam sizes for the radio images are shown 
at lower right (filled circle for the JVLA data, and open circles for 
the VLA data). The JVLA observations detected no radio emission 
within 28" of the cluster core. 

level of 1.58 /iJy beam -1 , within 9% of the theoretically- 
predicted thermal noise. The check source showed that 
coherence losses were < 1.5%. Thus, our 3a upper limit 
on the radio luminosity of Gl on 2011 June 26, assum- 
ing a flat radio spectrum extending up to the central 
frequency of 6.225 GHz, is < 2.0 x 10 31 ergs -1 . 

Given the results of UGH07, we can propose three pos- 
sible explanations for our highly-constraining radio non- 
detection. First, the emission could be extended (e.g. 
a pulsar wind nebula) and therefore resolved out in the 
higher-resolution JVLA observations. Alternatively, the 
radio source could be variable, by a factor of > 6, on 
a timescale of years. Finally, the detection reported by 
UGH07 could have been an artefact, arising either from 
the mix of VLA and JVLA antennas in use in 2006, or 
from the unfortunate coincidence of a noise spike with 
the known X-ray position. 

To test the first scenario, we tapered our JVLA data to 
the resolution of the C-configuration data of UGH07. No 
significant emission was detected within the 95% radio 
error circle, down to an rms noise level of 5.6 /iJy beam -1 , 
thereby ruling out extended emission. 

3.3. Reanalysis of archival VLA data 

To investigate radio variability, we reanalyzed the 
archival VLA observations of Gl (Table [2]). In the 
observations from 2006 and 2007, 5-6 of the original 
VLA antennas had been retrofitted with the new JVLA 
electronics. To account for the differing bandpasses 
of the retrofitted JVLA antennas, we applied baseline 
corrections within AIPS (jGreisenl l2QQ3h . after which we 
imported the data into CASA for further calibration 
and imaging, using the same procedures and calibrator 
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Table 2 

Radio observations of Gl. Where two different flux densities are quoted, the first is for the full data set, whereas the second excludes 

the newly-upgraded JVLA antennas. 



Program code 


Date 




Array config. 


Frequency 
(GHz) 


Bandwidth 
(MHz) 


Time on source 
(min) 


Source flux density 
(fj,Jy beam -1 ) 


Image noise 
(/iJy beam -1 ) 


AA 276 


2002 Sep 11 




B 


8.4 


100 


85 




47 


AG 730 


2006 Nov 24- 


26 


C 


8.4 


100 


845 


15/31 


5.5/7.2 


AK 634 


2007 Jan 13- 


14 


C 


4.9 


100 


214 


16/22 


AU 116 


2007 Jan 14- 


15 


c 


4.9 


100 


235 


67/63 


16/20 


SC 489 


2011 Jun 26 




A 


6.2 


2048 


446 




1.6 



sources as for the JVLA data (Section E2J). UGH07 
noted some issues with the JVLA antennas in the 8.4- 
GHz observations of 2006 November, which we found to 
be manifested as large (albeit smoothly-varying) phase 
corrections. As a conservative approach, we therefore 
made images both including and excluding data from 
the JVLA antennas from all observations from 2006 and 
2007. At 4.9 (but not at 8.4) GHz, there was sufficient 
emission in the field to apply phase self-calibration with 
a 10-minute solution interval. 

Only two data sets showed marginally-significant radio 
emission at the originally-reported position (Table [2]). In 
the 8.4-GHz observations of 2006 November, when ex- 
cluding the JVLA antennas we found a 4.3a source at 
31±7/iJy beam -1 , consistent with the results of UGH07. 
However, enhanced local noise in the source region (rms 
of 8.3 //Jy beam -1 ) reduces the true significance to 3.8a. 
The detection was not significant when the JVLA anten- 
nas were included. 

In the 4.9-GHz observations of 2007 January 14-15, 
there were no noticeable issues with the JVLA anten- 
nas. A 4.1a source was detected at a position consistent 
with the radio source reported by UGH07 (3.2a when 
excluding the JVLA antennas), although again the lo- 
cal rms appeared to be slightly enhanced. The data 
taken the previous day (January 13-14) showed no sig- 
nificant emission at the same position. Thus, either the 
radio emission is variable on timescales of 1 day, or the 
marginally-detected source is spurious. 

4. DISCUSSION 
4.1. Variability 

While our archival re-analysis confirms the measure- 
ment of UGH07, the enhanced local noise close to the 
source position reduces its significance. Coupled with 
the phase-center location and the slight po sitional off- 
set from the well-determined X-ray position (jKong et al.l 
l2010f ) , we cannot exclude the possibility that it is a noise 
spike. However, the probability of 4a noise spikes co- 
inciding in two independent images (Fig. [1]) makes this 
interpretation unlikely. 

Should these two marginal detections be real, they 
would imply a variation in the radio emission by a fac- 
tor of > 6 between 2006 and 2011, and variations on 
timescales as short as 1 day in 2007 January. Radio emis- 
sion at this level (~ 10 32 ergs -1 ) could correspond to 
sub-Eddington flaring from an IMBH, or to a major out- 
burst of a black hole LMXB (an outbursting neutron star 
would not reach the observed radio luminosity). Stellar- 
mass black holes are known t o exhibit short-timescale , 
sub-Eddington radio flares (e.g. iMiller- Jones et al.l l2008). 




1.638 log L fl - 1.136 log L x - 6.863 

Figure 2. Our best-fitting regression for the FP (Equation [T} as 
a function of black hole mass, showing the limits implied by our 
simultaneous radio and X-ray observations of Gl. 

Should the radio emission be an analogous flare from 
an IMBH, then scaling the variability timescale with ac- 
cretor mass would imply radio variations over weeks to 
months rather than days, making this scenario unlikely. 

The more likely explanation is an LMXB outburst. 
Given the observed correlation between outburst du- 
ration, peak luminosit y , and o rbital period in LMXBs 
(|Portegies Zwart et al.l |2004b), the radio luminosity 
would imply an extended period (several months) above 
10 36 ergs _1 . Bright radio flares occurring over a period 
of a few months and varying on times cales of days have 
been seen in black hole LMXBs (e.g. McClint ock et all 
2009). The lack of X -ray coverage between 2003 and 
2008 implies that the corresponding X-ray flare would 
not have been detected. 

4.2. Mass constraints from the Fundamental Plane 

Although black hole masses can in principle be 
constrained using the FP relationship (Mer loni et al.l 
2003), the error bars on the best-fitting regression co- 
efficients are relatively large. Furthermore, the re- 
gression is sample-dependent, with the intrinsic scat- 
ter about the FP being minimi zed when using onl y 
strongly sub-Eddington objects (Kording et al. 2006). 
While most published FP rela tions ( Merloni et al.l 120031 : 
iKording et al1l200fl : iPlotkin et aTll2012ft have used radio 
or X-ray luminosity as the dependent variable, the intrin- 
sic scatter about the plane implies that it is not valid to 
simply invert the best-fitting regression coefficients to de- 
termine a mass constr aint from the me asured radio and 
X-ray luminosities (Plotkin et al. 2012). The sole exist- 
ing mass regression ([Gultekin et al.l 120091 ) included only 
supermassive black holes (M > 1O 6 M ), so has limited 
predictive power for lower-mass systems. 

To place a valid constraint on the mass of the X-ray 



No radio emission from an IMBH in Gl 



5 



-0.5 
logM Bffjme 



Figure 3. Histogram of the residuals from our best-fitting FP re- 
lation (Equation^), together with a Gaussian of the same mean and 
standard deviation (0.44 dex). The scatter about the best-fitting 
regression implies that black hole masses cannot be predicted to 
better than a factor of ~ 3 using the FP, and that FP mass pre- 
dictions can be in error by over an order of magnitude. 

source in Gl, we therefore perf ormed a new regres sion of 
the sub-Eddington sample of iPlotkin et al.l (|2012f h with 
mass as the dependent variable. The inclusion of BL 
Lacs in the sample (having controlled for synchrotron 
cooling effects) makes the implicit assumptions that the 
radio emission has a flat or inverted spectrum, and that 
the X-rays are predominantly optically-thin synchrotron 
emission. However, it improves the dynamic range in 
black hole mass, and any systematics introduced by rela- 
tivists beaming are of secondary importance when com- 
pared with th e intrinsic scatter about the FP. Following 
IPlotkin et al.l (j2012[ ). we used the Bayesian technique of 
iKellvl (|2QQ7| ) to handle the coupled uncertainties in radio 
and X-ray luminosities. We find a best-fitting regression 

log Mbh = (1.638 ±0.070) log L R 

- (1.136 ± 0.077) logLx - (6.863 ± 0.790), 

(1) 

where Lr and Lx are the radio and 0.5-10 keV X-ray 
luminosities, respectively, in ergs -1 . Our simultane- 
ous measurements then predict a black hole mass of 
< 1.7 x 1O 3 M in Gl (Figure EJ). Comparing measured 
black hole masses with those predicted by the regres- 
sion shows that the la uncertainty on the mass esti- 
mates is 0.44 dex (Figure EJ). This is significantly higher 
than the 0.12 dex used by previous authors to determine 
the uncertainties on FP-determined IMBH masses (e.g. 
iCseh et al.ll2010UKirsten fc Vlemmings!l2012[ ). but, being 
empirically-determined from the data, should be a more 
accurate estimate of the true uncertainty. Thus, we esti- 
mate the 95% confidence upper limit on the mass of the 
X-ray source in Gl, if a black hole, to be < 9.7 x 1O 3 M . 

4.3. Does the X-ray emission arise from an IMBH? 

iKong et "all (j2010f ) concluded that the only possible 
explanations for the X-ray emission from Gl were an 
IMBH accreting from a companion star, or a typica l 
LMXB. As discussed by iPoolev fc Rappaportl (j20Q6h . 
seven persis tent LMXBs are known in Galactic globu- 
lar clusters (|Verbunt fc Lewinl l2006f ) at X-ray luminosi- 
ties > 10 36 ergs _1 . Since Gl is a particularly mas- 



sive, dense cluster (jMevlan et al. 2001), it likely contains 
many LMXBs. The stability of the X-ray luminosity be- 
tween 2001 and 2011, coupled with the positional coinci- 
dence between the Chandra detections of 2008 and 2011, 
makes it unlikely (albeit not impossible) that some frac- 
tion of the X-ray emission arises from multiple LMXBs 
in the cluster core. However, should LMXBs contribute 
some or all of the observed X-ray luminosity, then our up- 
per limit on the mass of an IMBH in Gl would increase, 
via Equation [TJ 

Given the observed X-ray count rate, there is no way 
to distinguish the X-ray spectrum of an IMBH from an 
LMXB with current instruments. Combining the lack of 
strong radio emission from Gl with the fact that LMXBs 
of the observed luminosity and spectral shape are rela- 
tively common in globular clusters like Gl, an LMXB 
origin for the X-ray emission is the most probable inter- 
pretation. However, given the positional offset from the 
marginal radio detections of 2006-2007, the radio and 
X-ray emission likely arise from different LMXBs. 

4.4. Intermediate-mass black holes in globular clusters 

Although IMB Hs are predicted to exis t in globu- 
lar clusters (e.g iMiller fc Hamilton! l2002f ). many at 
tempts to de tect central dark ma s ses via stella r kine- 



matics (e.g. iGerssen et all 120031: iNovola et all [2008) 
have been controversial ( see Baumgardt et al.l [2003a; 



Ivan der Marel fc Anderson! 120101 . respectively). While 
deep radio observations have the potential to discrim- 
inate between an IMB H and a collection of dark rem- 
nants (|Maccaronell2004l ) , our result implies that no radio 
observations to date have shown conclusive evidence for 
an IM BH in a globular cluster (Mac carone fc Servillatl 
2008). The most stringent radio limits available imply 
that either IMBHs are rare in globular cl usters, or that 
they are extremely inefficient accretors (S tr ader et al.l 

[ml). 

5. CONCLUSIONS 

Deep radio continuum observations with the JVLA de- 
tected no radio emission within 28" of the X-ray source 
at the center of Gl, to a 3a level of 4.7 /iJy beam -1 . The 
0.5-10 keV X-ray luminosity measured simultaneously 
with Chandra was 1.74^-j 3 x 10 36 (d/750 kpc) 2 ergs -1 . 
Using these measurements together with a new FP re- 
gression, we constrain the mass of the X-ray source in Gl, 
if a black hole, to be < 9.7 x 1O 3 M , at 95% confidence. 
Our 3<j radio upper limit is a factor of 6 deeper than 
the previously-reported VLA detection (UGH07), sug- 
gesting either that the previously-detected source was an 
artefact, or, more likely, that the radio emission is time- 
variable, and arises from a black hole LMXB in outburst. 

While the FP correlation works well for a large sam- 
ple of sources, it is a statistical tool, and as Fig. [3] 
shows, can be in error by over an order of magnitude 
for individual sources. Therefore, despite our radio 
non-detection, we cannot definitively exclude the X-ray 
source in Gl being an IMBH. T hus, the 4<r dynamical 
detection (jGebhardt et al.ll2005f ) still makes it a viable 
globular cluster IMBH candidate. Our current radio lim- 
its are unlikely to be surpassed in any reasonable inte- 
gration time with the current generation of telescopes, so 
unless the possible variability can be confirmed, a defini- 
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tive identification of the X-ray source will remain the 
preserve of future instruments. 
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